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We study the transformation of multipolar decompositions of gravitational radiation under ro- 
tations and boosts. Rotations to the remnant black hole's frame simplify the waveforms from 
the merger of generic spinning black hole binaries. Boosts may be important to get an accurate 
gravitational- wave phasing, especially for configurations leading to large recoil velocities of the rem- 
nant. As a test of the formalism we revisit the classic problem of point particles falling into a 
. Schwarzschild black hole. Then we highlight by specific examples the importance of choosing the 

' right frame in numerical simulations of unequal-mass, spinning binary black-hole mergers. 
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I. INTRODUCTION 



^ ' Multipolar expansions are a fundamental tool in the study of gravitational radiation. Such expansions can be 
qh| performed using different sets of basis functions (see [J for a review). It is common practice in numerical relativity 
^ simulations to decompose the Weyl scalar ^4 (a quantity related to outgoing gravitational radiation) as a sum over the 
j^jT)^ Newman-Penrose spin-weighted spherical harmonics (SWSHs) of spin-weight s, commonly denoted by sYim- SWSHs 
reduce to ordinary spherical harmonics for s — 0, and |s| = 2 is the spin- weight of interest for gravitational radiation 
ii. 

Compact binaries are one of the most promising candidate sources of detectable gravitational waves for Earth- 
' based and space-based interferometers. The post-Newtonian approximation, an analytical expansion of Einstein's 
^— I ' equations valid for orbital velocities F/c <C 1, is usually employed to describe the early inspiral of compact binaries. 
■ The post-Newtonian expansion becomes increasingly inaccurate late in the inspiral, when orbital velocities become 
' large. At these late stages, time evolutions of compact binaries in numerical relativity are necessary. High-accuracy 
. numerical evolutions are difficult for various reasons, including the approximate nature of the initial data, errors 
' arising from discretization and boundary effects. Post-Newtonian expansions and numerical calculations must be 
00 , cross-validated by systematic comparisons of the resulting waveforms, in order to determine a common region of 
' validity and to build reliable phenomenological templates [4, 5, 6, 7,, ^, .9., JJI, .11, lISi] ■ Since numerical waveforms are 
J> ] usually decomposed in SWSHs, significant efforts have recently been devoted to study multipolar SWSH expansions 
both within post-Newtonian theory [l^, [l3| and in connection with numerical simulations [3, [IB, [13, [H, Il9l |. 

This paper deals with the transformation properties of SWSH decompositions of gravitational radiation under 



■ rotations and boosts. We will mostly focus on the radiation emitted during the inspiral, merger and ringdown of 
■ ■ ' compact binaries, but our main results are general enough that they should be useful in more general situations. 

The transformation properties of gravitational waves have been studied as part of the general problem of spacetime 
transformations in the sixties 0, [20 . [2ll . [2^ . It was realized that the asymptotic symmetries of asymptotically 
flat spacetimes form an infinite-dimensional group, now known as the Bondi-Metzner-Sachs (BMS) group, and that 
Lorentz transformations are a subgroup of the BMS group. Later on, the transformations of SWSHs under the BMS 
group have been investigated in a classic paper by Goldberg and collaborators [^ . 

Unfortunately these pioneering works do not give an explicit prescription to transform gravitational waves under 
rotations and boosts, and to handle their multipolar decomposition. Our main purpose is to fill this gap and to 
provide examples of the formalism at work, by considering specific problems in numerical relativity and black-hole 
perturbation theory. We only consider the emission as seen by an observer far away from the source, and we assume 
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the spacetime to be asymptotically flat. As discussed in the next section, the spacetime can then be described by a 
set of coordinates {u,r,9,(j)), where 9 and (f) are coordinates on the sphere, m is a retarded time parameter and the 
limit r — !■ cx) (with u finite) represents future null infinity. In this limit r is much larger than u and any translation 
Sx of the source: |i5a;|/r <C 1, \u\/r <C 1. Furthermore, when we consider boosts we will always neglect quantities of 
order (V/c)'^, where V is the boost velocity of the source. This approximation is well justified: according to recent 
results from numerical relativity, black hole mergers produce recoil velocities V < 4000 km/s, i.e. V/c < 0.013 

In the remainder of this introduction we list some of the main motivations for our work, and provide a brief summary 
of our most useful results. 



A. Motivations 

The multipolar structure of gravitational radiation is an important element in the problem of detection and pa- 
rameter estimation from compact binary inspirals. The energy distribution determines which {l,m) components of 
the radiation are most easily detectable. For a binary inspiral, it has long been believed that the I = m = 2 mode 
should dominate the energy emission [2^. Numerical simulations confirm that the I — m = 2 mode carries ~ 98% 
of the radiation for equal-mass, non-spinning binaries (isj . However, higher modes can carry more than ^ 10% of 
the energy for unequal-mass binaries [IGj and/or when the black holes are spinning [3,[l^. Now, recall that the 
gravitational-wave energy scales with the square of the waves' amplitude. Since SWSHs are normalized over the 
sphere, this means that the amplitude ratio of the I — m — 3 and I = m — 2 modes is typically « 1/3 for generic 
binaries when different multipolar components are averaged over the sky. Depending on the relative direction of the 
source and the detector, the enhancement of multipoles with I > 2 can be even larger (see [2^ for a discussion in 
the context of ringdown detection). For generic binary inspirals, a good phase accuracy on the I — m — 2 mode 
is necessary but not sufficient to produce accurate templates: higher multipolar contributions must be taken into 
account to carefully track the phasing of the waves over many inspiral cycles. One may expect that fewer multipolar 
components should be sufficient in the early inspiral phase, but this expectation is not confirmed by calculations. For 
example, by comparing multipolar expansions of the (numerically computed) adiabatic energy flux to post-Newtonian 
expansions of the flux in the context of extreme mass-ratio inspirals, Ref. 12] showed that the number of multipoles 
required to achieve a given post-Newtonian accuracy actually increases in the Newtonian regime (i.e. for low orbital 
frequencies) . 

The transformation of multipolar decompositions of gravitational radiation under rotations deserves a detailed 
study for many reasons, some of which we list below. 

1) Provide translation rules between numerical codes using different conventions. A pressing motivation for our 
work is of practical nature. As more accurate and complex numerical waveforms are produced, it is vital to establish 
a set of conventions between different numerical groups, and rules to translate between different conventions. Binary 
simulations often differ in subtle details of the numerical settings: for example, some groups initially place the black 
holes on the y-axis, while other groups place them on the a:-axis. This is obviously a matter of convention, but 
conventions are important when setting up data formats to be used in LIGO and Virgo data analysis (see e.g. (30j 
for an attempt in this direction). Similar issues should also be taken into account in collaborative efforts to inject 
numerical waveforms in the LIGO data stream, such as the Ninja project [sij. One purpose of this paper is to provide 
a "dictionary" to translate between different conventions by relating multipolar expansions in different frames. 

2) Identify the optimal frame to study the merger /ringdown phase in generic spinning binary simulations. Given a 
generic configuration for a spinning, inspiralling binary, at present there is no reliable analytical method to predict 
(given the initial data) the orientation of the remnant's spin. We will refer to the frame where the initial orbital 
angular momentum points in the z-direction as the simulation frame, and to the frame where the remnant black 
hole's spin points in the ^-direction as the remnant frame. 

In the simple (and most-studied) case of non-spinning binaries, a distinction between the simulation frame and the 
remnant frame is not necessary. The black holes usually orbit on the x — y plane, the initial orbital angular momentum 
is in the ^-direction, and angular momentum is radiated only in the z-direction. Then the spin of the final black 
hole will also point in the z-direction, and its dimensionless magnitude j — J/M'^ can be computed by evaluating the 
z-component of the radiated angular momentum and using balance arguments. Since the remnant's spin points in 
the z-direction, a fit of the ringdown frequencies for the dominant SWSH component of the waves in the simulation 
frame can then be used to verify the accuracy of wave extraction methods 

For generic inspirals, however, the system will radiate angular momentum in the x- and y-directions as well. The 
spin of the final black hole will not be aligned with the z-axis of the simulation frame. If we (incorrectly) perform the 
SWSH decomposition in the simulation frame, rather than in the remnant frame, the expansion will mix quasinormal 
modes with the given value of I and different values of m. For this reason, using the wrong frame produces errors in 
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the quasinormal frequencies of the final black hole and in ringdown-based estimates of the value of the final angular 
momentum. Note that modes with different values of I are decoupled, since SWSHs with fixed I and s form an 
irreducible representation of the rotation group. In Section IIVI we will illustrate this problem by specific examples, 
and we will show that rotations can be used to simplify the waveforms: as viewed from the normal to the final black 
hole's rotation plane, the merger-ringdown radiation in the remnant frame is circularly polarized. 

3) Measure the spin direction of the remnant through merger /ringdown waves. The measurement of a single complex 
quasinormal frequency with Earth-based and space-based interferometers [s^ allows a determination of the mass 
and spin magnitude of the hole resulting from a merger (the I = m = 2 fundamental quasinormal mode is again 
expected to be dominant). For Kerr black holes all quasinormal frequencies are uniquely determined by the mass 
and spin of the hole, and the measurement of the next subdominant mode (typically the fundamental mode with 
Z = m = 3) will provide a test of the Kerr nature of the merger remnant. Since different multipolar components have 
different angular dependence, a multi-mode detection could even be used to measure the spin direction of the remnant. 
A detailed analysis of this problem requires the knowledge of the transformation properties of the multipoles under 
rotations and boosts, and it is an interesting topic for future work. 

As argued earlier, boost velocities in binary mergers are typically small [V/c < 0.013). However the study of boosts 
is of more than academic interest. In numerical simulations, gravitational radiation should ideally be extracted on 
a sphere at null infinity. In practice the extraction sphere is usually located at sufficiently large, but finite, distance 
from the source. In this paper we point out that the distortion of the radiation pattern induced by a large kick of 
the remnant can induce a mode mixing in the signal. To leading order the amplitude of this mixing is proportional 
to the kick velocity. The mixing effect will typically be small iV/c ~ 10^^ for the largest predicted kick velocities), 
but it may still affect the phasing error of the inspiral waves over long inspiral times. Conversely, the magnitude of 
the boost-induced mode mixing may set a threshold on the number of high-order post-Newtonian terms we can (or 
should) compute to get an accurate phasing. Finally, one should not exclude a priori the possibility that the large 
recession velocities due to Hubble flow may produce observable effects on gravitational waveforms. 

Boosts are particularly relevant for the extraction of gravitational waves from the collision of highly relativistic black 
holes with relatively small impact parameter. This process is of interest as an extreme testing ground of strong-field 
gravity and in the context of mini black hole production at the Large Hadron Collider, and it is currently under 
investigation [ssj . 

B. Main results 

For the sake of clarity, and in view of possible generalizations of our results, in the body of this paper we give a 
detailed account of our calculations. To help the reader, here we provide a short list of results that we consider most 
useful in phenomenological applications: 

• Eq. (|3.23p relates the Lorentz transformation of spin-2 SWSHs to the transformation of the Weyl scalar ^'4. 

• Eq. p.28p can be used to obtain spin-2 SWSHs in stereographic coordinates for general values of I, and Eq. (|3.30p 
gives closed- form expressions valid for I = 2 and Z = 3. 

• For pure rotations different V s do not couple because SWSHs with given (s , I) form an irreducible represen- 
tation of the rotation group, and we can consider the transformation laws separately for each value of I. The 
transformation of the harmonics under general rotations, in terms of a constant rotation matrix, is given in 
Eq. 

• Eq. p.53p is one of our main results: it provides an expression for the general rotation matrix, valid for all 
values of I. Eqs. (|3.46p and (|3.54p list components of the rotation matrix for / = 2 and / = 3, respectively. 
Eqs. (|3.49p and (|3.56p further specialize these results to frames related by a 7r/2 rotation (i.e. 9 ^ 9' = 9 — 7r/2) 
The transformation under rotations of the Weyl scalar ^'4 and of its complex conjugate has a particularly simple 
form, which is given in Eqs. (j3.59b[) . (|3.60b[) . 

• Eq. (j3.58p provides the I — 2 transformation matrix under non-relativistic boosts, and Eq. p.65p shows in a 
special case how boosts produce mode-mixing. As far as we know, this fact has never been pointed out in the 
existing literature. 

• In Section [IV Al we show that our transformation matrices are consistent with two "classic" calculations. In the 
first calculation we consider the equal-mass head-on collision of black holes as seen in two frames related by a 
7r/2 rotation (i.e. 9 ^ 9' = 9 — 7r/2). Then we look at particles falling radially into Schwarzschild black holes 
along the z-axis or along radial, equatorial geodesies. We show that the SWSH components of the stress-energy 
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tensor in these two frames are related by proportionality constants which are consistent with our transformation 
laws. 

• In Section IIVBI we demonstrate that our work can be useful to study the generic merger of unequal-mass, 
spinning black holes. We suggest a simple recipe to switch from the simulation frame to the remnant frame. 
Then we show that the merger waveform does indeed look simpler in the remnant frame, because in this frame: 
(1) ringdown radiation is circularly polarized, and (2) ringdown frequencies are consistent with the final black 
hole spin predicted from balance arguments. These properties do not apply to the simulation frame. 

The paper is organized as follows. In Section[n]we introduce some basic notions on the BMS group. Section|TTT]is the 
core of this paper, where we derive our main results on the transformation properties of SWSHs and of gravitational 
radiation under rotations and boosts. Section IIVI applies the formalism to specific problems in numerical relativity 
and perturbation theory, and Section |V] discusses some directions for future research. The Appendices contain various 
technical topics: the expression of non-relativistic boosts in the context of BMS transformations (Appendix IX|: the 
description of rotations and boosts in the stereographic representation (Appendix [B| ; a discussion of our numerical 
method to compute the radiated angular momentum for generic, spinning black-hole binaries (Appendix [C| ; and a 
brief summary of perturbation theory calculations of particles falling into black holes (Appendix iPj) . Unless mentioned 
otherwise, throughout the paper we use geometrical units {G = c = 1). 



II. THE BONDI-METZNER-SACHS GROUP 

To begin with, in this section we briefly recall the main steps of the BMS-group construction. We refer to [3, [111, 
Ib^ . [s^ . 36] for a more detailed and rigorous account. Any asymptotically flat spacetime can be described by standard 
coordinates {u,r,9,4>) defined as follows [13]: 

• M is a retarded time parameter, defined by requiring that the hypcrsurfaccs u ~ constant are null: g^'^u^fj.u^i, — 0. 
These hypersurfaces are the wavefronts of gravitational waves emitted by the source. 

• r is an affine parameter of the null geodesies tangent to = g^'^u^u {l^ is also a vector of the Newman-Penrose 
tetrad) . 

• 9 and cj) are polar coordinates on the manifolds (diffeomorphic to the 2-sphere) r — constant, u = constant. 
Geodesies with constant (it, 9, (j)) are called "rays," and they are a congruence of geodesies, generators of the 
hypersurfaces u = constant. The limit r —^ oo with finite u corresponds to future null infinity, where the 
gravitational radiation emitted by the source is detected. 

• Far away from the source one can define a time coordinate t = u + r. 

The spacetime is asymptotically flat if the metric tends to the flat space metric sufficiently fast as r ^ oo (see 
(2l| for further details). This construction allows one to define asymptotic symmetries and then, quoting Sachs, to 
"separate the kinematics of space-time from the dynamics of the gravitational field at least at spatial infinity" [2l[ . It 
is worth noting that such a separation can only be done in asymptotically fiat spacetimes, by looking at asymptotic 
symmetries; otherwise, "in general relativity such a separation is usually impossible" |21i] . 

Asymptotic symmetries are the transformations preserving the asymptotically fiat form of the metric, i.e. pre- 
serving the above properties. These symmetries form a group: the BMS group. The vector fields generating such 
transformations satisfy an equation which is weaker than the Killing equation, and coincides with the Killing equation 
as r — > oo. The BMS group is defined by the following transformations of the (w, 9, (j)) coordinates: 

lim 9' = 9' [9,(1)), 

r— i-oo 

lim (/)' = <j)'{9,<j)), 

Mm u' ^ K{9,(t))[u~ a{9,(t))], (2.1) 

r — *oo 

where 9' {9, cf)), 4)'{9, (j)) and K{9, (f) are related by a conformal transformation, i.e. 

d9'^ + sin^ 9'd(l)'^ = {d9^ + sin^ 9d(t>'^) . (2.2) 
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We follow the convention of Newman and Penrose [3], but note that Ref. Q uses a different convention: their K is our 
K~^. We can consider Eqs. (|2.ip as transformations defined on a three-dimensional manifold described by {u,9,(p), 
i.e. on the future null hypersurface: 



u' = K{0,(l))[u~ a{9,(t>)]. (2.3) 

The transformations apply to asymptotic quantities defined at future null infinity, which are obtained by expanding 
tensors (or Newman- Penrose scalars) in powers of 1 /r (STj . This approach is particularly meaningful in the study of 
gravitational-wave emission, since the detector's location is usually a good approximation to future null infinity. 

The BMS group is the semidirect sum of the Lorentz group with the group of supertranslations, which are a 
generalization of translations. The supertranslations are the transformations (|2.3p with 9' = 9 , (j)' = (j) and K{9, 0) = 1. 
A useful illustration of their meaning is given in Ref. [34]. Translations are special cases of supertranslations, with 
a{9,(j)) = Co + tiTii. Here « = 1,2,3 and n,; = Xi/r = (sin cos sin sin 0, cos 0) are director cosines. 

Lorentz transformations are transformations of the form (|2.3p with a{9, 4>) = 0: 

9' = 9'{9,<l>), 
cf>' = <ly'{9,4,), 

v! = K{9,(f>)u. (2.4) 

Strictly speaking, there is no unique way to identify the Lorentz group as a subgroup of the BMS group 0, Is^ . 
Indeed, if {L} are the transformations (|2.4p and T{a) is a supertranslation with fixed a, then the transformations 
{T{a) LT~^(a)} form another group isomorphic to the Lorentz group. As far as we know, there is no rigorous proof 
that our choice (j2.4p of the Lorentz subgroup is more "natural" than others, but there are arguments supporting 
this choice. In particular, Ref. [1] suggests that with further, reasonable physical assumptions, the Lorentz group 
embedded in the BMS group is restricted to {T{a) LT~^{a)} with a = eo + CiUi (ordinary translations); the residual 
freedom, then, corresponds to an arbitrariness in the choice of the origin. Therefore it is not restrictive to assume 
a = 0, i.e. to assume that the Lorentz transformations are given by Eq. (|2.4p . Ordinary rotations are special cases of 
Lorentz transformations with K{9, (f) — 1. As discussed in Appendix [XI a boost with (non-relativistic) three- velocity 
V = (Fi, V2, Va) has K{9, </>) = ! + V,n,. 

To study Lorentz transformations, and then conformal transformations on the 2-sphere, it is useful to introduce 
the stereographic representation, in which the coordinates {9 , 0) on the 2-sphere are mapped on the complex plane 
by introducing a complex variable 

C = cot ^ e"^ . (2.5) 
Lorentz transformations correspond to mappings C — > C' of the form 

where a,b,c,d are complex numbers and we can assume without loss of generality that ad — be = 1. The interested 
reader can find a more detailed discussion in Appendix [Bl 



III. BMS TRANSFORMATIONS OF GRAVITATIONAL RADIATION 



As we discussed above, = g^'^u^i, is a null vector, and r is an affine parameter of the null geodesic tangent to 
l^^ . Following 0, [33, [3§|, we normalize r by requiring that = dx^^/dr; furthermore, we take ^'^ to be a vector 
of the Newman-Penrose tetrad. It is important to stress that this is a particular choice of the tetrad [13]: we are 
fixing the tetrad boost freedom AZ^, n'' X~^n'^ (and then the normalization of ^4, since for a tetrad boost 

^'4 A^^^4). Ref. [131 calls this choice the "special coordinate system and associated tetrad", and Ref. ^] calls 
it the "Bondi tetrad." Most of the literature on the asymptotic behavior of Newman-Penrose quantities and on 
BMS transformations assumes this tetrad choice. The same normalization has been used by Teukolsky to study 
perturbations of Kerr black holes [4l[ . As discussed in |4^ , recent works dealing with numerical evolutions of the 
Einstein equations use a tetrad choice such that la = u^al\f^ [l^, [IB, IS IS] ■ We do not follow the latter convention 
because Eq. (|3.3p below, as derived in [sj, requires la — u^a- 
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A. Weyl scalar and shear 

Gravitational radiation is described by the Weyl scalar '^4. The asymptotic limit of ^'4 at future null infinity is 
closely related to the asymptotic limit of the shear a, which is defined in terms of the Newman-Penrose tetrad as 

cr(u, r, 6, (j)) — If^-^ijm^my . (3-1) 

Indeed, assuming the special coordinate system and the associated tetrad defined above, the asymptotic expansion of 
^'4 and g is [37II 

r 

air,u,9,^) = + o{r-') , (3.2) 

where 

and an overbar denotes complex conjugation. To simplify calculations, in the remainder of this paper we will usually 
consider SWSHs of spin weight s = 2 and expand the complex-conjugate Weyl scalar as: 

^0 ^ g-i^n ^ 2 V'im 2Yl^ 0) , (3.4) 
Im 

where we assume a gravitational- wave solution with given frequency lo. In numerical relativity simulations the Weyl 
scalar is usually expanded in SWSHs of spin- weight s = — 2. The relevant expansion can easily be obtained by recalling 
that SWSHs satisfy the property 

sYiUO, <P) = (-1)"+^ ^sYi _„(0, 0) , (3.5) 

which implies that 

^0 ^ gi^u ^ -2Ylrn{e, <f) (3.6) 

Irn 

with 

_2Vim= (-l)'"2V5i-m. (3.7) 

The first derivative of the shear with respect to retarded time, da^ /du, is called the news function. Its squared 
modulus represents the energy density of the gravitational wave: the total energy emitted by the source (corresponding 
to its mass loss) is 







An j 


du 



2 



dudVL . (3.8) 



To understand the relation of ct'^ with the gravitational-wave amplitude, let us consider the asymptotic metric in the 
transverse-traceless (TT) gauge; following the notation of Teukolsky [4l|] , 

*4 = ^(/i+-i/ix) = ^, (3.9) 

where we have defined the complex wave amplitude h in terms of the two polarization modes and as h = h-\-~\h^ . 
Assuming again a gravitational-wave solution with given frequency, 

^ (3.10) 

r 



we have 

2du 



*4 = 7TTr^. (3.11) 



7 



From Eq. (|3?3)) it follows that 



rh = e"^"/i 



therefore, ct° corresponds (modulo proportionality constants) to the wave amplitude h. 
From Eq. (|3.8p we have 



1 



2^2 



167r 



■(|/i+P + |/i> 



which is the well-known expression for the energy carried by a gravitational wave |44l | . Notice that 



d^E _ r^\^l\ 



dtdil Attlu^ 



(3.12) 



(3.13) 



(3.14) 



which is consistent with Ref. [iH. In [HE Em, where la = u^a/V^, the Weyl scalar is ^'4 = {h+ — ihx)- With 
their normalization, d^E/dtdVl — r^j^'H/lGTra;^. The interpretation of the asymptotic shear as the amplitude of the 
wave in the TT-gauge is also discussed in Ref. fisj . 

Under a BMS transformation the asymptotic shear transforms as 0, [22] 



where the operator ^, acting on a quantity rj of spin weight s, is defined by 

(sin 0)^*77 



(3.15) 



(3.16) 



and the angle x represents the rotation of the vectors m'', to^ of the Newman-Penrose tetrad. For Lorentz transfor- 
mations a = 0, and Eq. p.lSp reduces to 



a°'(ii, 



(3.17) 



Incidentally, this relation also holds for ordinary translations, since p^a — in this case. The total radiated energy 
(|3.8p transforms as 







An J 


du' 



du'dn' ^ — I K 
47r 







1 K-^ 


du 



Kdu K^dn = — I K 
in 









du 



dudQ , 



(3.18) 



with a shift factor K~^, as expected. This relation is quite interesting: in the case of a boost, it describes the beaming 
of the emitted gravitational radiation. 

The transformation p.l7p is not yet complete, since cr^' is still expressed in the unprimed coordinates {u, 9, (p). The 
complete transformation involves the relation between primed and unprimed coordinates: 



a°'(w', e\ 0') = X-ie2'x(«(«')''^('^'»(T0(u(H'), 0(0')X0')) • 
For a gravitational wave with frequency oj we have (j'^{u, 6, (p) — e~"^"(7°(6', (p), and since u' = Ku, 

(tO'(u', e', (j)') = i^-ie2'>^(''(''')'"^(^'»e-'-^"'aO(0(0'), -^(0')) • 



(3.19) 



(3.20) 



We see that the frequency of the gravitational wave is blueshifted or redshifted: uj ^ uj' — uj / K . 

The transformation of the asymptotic Weyl scalar ^1 (we consider the complex conjugate ^"4, with spin 2, instead 
of 4*4, with spin —2, to simplify calculations), related to cr° by Eq. (|3.3p . immediately follows: 



^':iu',e',(b')=K-'e 



^liuiu'),9{e'),cP{<l>')). 



We now expand the complex-conjugate Weyl scalar in SWSHs of spin weight s = 2, as in Eq. 
transformation of the harmonics 



(3.21) 
Given the 



2W(^'(e,'^),0'(0,</>)) = ^ A ''in' dm. 



(3.22) 
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the Lorentz transformed Weyl scalar is 

I'm' 

with 

2'4''l'm' = E 2^'™^r™;i'm' • (3-24) 

In the case of rotations K = 1 and Dim-Vm' is block diagonal in I, because spherical harmonics with fixed I and s 
form irreducible representations of the rotation group. Then the transformation of the Weyl scalar reduces to 

Im' 

I 

2^L' = E 2i>imD%-}, (3.25b) 

m——l 

where we have defined 

D^mra' ^ A™;;™' • (3.26) 

Therefore, the transformation of the multipolar decomposition of the Weyl scalar under rotations is given by the 
(inverse) transformation of s = 2 SWSHs, modulo an overall phase factor. 

We now turn to a detailed treatment of the transformation of SWSHs under BMS transformations. We start with 
the general case, then specialize to rotations and non-relativistic boosts. Finally we show that there is a simple 
relation between the transformations of the Weyl scalar components in a SWSH expansion, and the transformations 
of SWSHs themselves. 



B. Transformations of spin- weighted spherical harmonics 

The BMS transformations of SWSHs are discussed in Ref. Given s and / > |s|, choose arbitrarily an integer 
L > I. Then SWSHs can be written as 

^^™=E E (3.27) 
m— 712—0 

where s^^^^"^ are constants whose explicit values can be found in Ref. and sZ^^^^ are properly defined functions 
of Q and This equation can be rewritten in a simpler form: 



lYim = (-1)' 



21 + 1 {I + my.il - my. 
(/ + 2)!(;-2)! 



E(-i)" 



1-2 
n 



1 + 2 

I — m — n 



where max(0, —2 — m) < n < I ~ max(m, 2) and we have defined 



(3.28) 



(3.29) 



This expression is consistent with Eq. (4.13) of Ref. i3;] if we correct for an erroneous factor of (—1)™ in their Eq. (2.12). 
For ; = 2 and / = 3 the SWSHs reduce to 



25^2m — (^1) 



1 



TT ^(2 + TO)!(2-m)! 



2 J^3m 



1 - 5,n.z jS-my _ 1 - S„i-3 / (3 + to)! 
(2-m)!V (3 + m)! " (2 + m)! V (3 - m)! 



(3.30a) 
(3.30b) 
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where we have defined 



2+m 

2^oV2 = (^:^C'-™= (^sin^j [cos"-j e'™*e-2'^ (-2 < m < 2) , (3.31a) 
^ = ^^-—C^-™ = (^sin - j (^cos - j e^-^e-^'-* , (-3 < m < 2) , (3.31b) 



1 


(1 






1 


(1 






cc 


(1 





lYi^n = , , ( - cot dde + —2- - — (ae - cot 0) ) F,™ , (3.32) 



24-m / /J \ 4~m 

^ 2^0^^^., = ^^Jil^C'-™ = ( Sin ^ j (^cos - j e'™*e-2"^ , (-2 < m < 3) . (3.31c) 

The factor e^^"^ is due to different conventions in defining SWSHs, corresponding to different choices of the Newman- 
Penrosc tetrad (see footnote 6 of [s*]). We have checked that Eqs. p.28p and p.30p are consistent with the standard 
definition of SWSHs in terms of ordinary spherical harmonics: 

2m 

^(/- l);(? + l)(/ + 2) V" " ' sin^e sm 

where we have foUowed the conventions of j3i] , adopted also in recent numerical simulations [UK 113, El. 

The only 

difference between Eq. p.32p and our expressions is a global factor of e~^"^, which does not depend on to and is 
therefore irrelevant. The s = —2 harmonics can be obtained by complex conjugation, using Eq. p.Sp . 

Let us focus for the time being on i = 2. In stereographic coordinates (see Appendix IB|) a Lorentz transformation 
maps C ^ C' as in Eq. (|2.6p . and 

(aC + c)^-"' (6C + rfT+" ^K'e^^ ^ F^-^Z™ , (-2 < to < 2) , (3.33) 

where the constants Fm'm are given by the expansion coefficients of polynomials of the form 

2 

(aC + 5)2-"' (6C + rf)' = Prn'm'e-'^ , ("2 < to' < 2) . (3.34) 

m=-2 

Replacing Eqs. (|3.33p and (|3.34p in Eq. p.30ap we find the transformation of spin-2 SWSHs with I = 2 and its inverse: 

2 

2ix(e.0) /l(2) 



if2(0,0)e2ix(M) ^ Al]^2Y2m(e,<f), (3.35) 

m=-2 
2 

>2™(0,0) = ;f-2e-2ix(fl',0') ^ A^^^^,\Y2m'{e'A'), (3-36) 



m' = -2 

where 



,(2) _ f ..m+m' I (2 + to)!(2-to)! 
^m'm-l y(2 + TO')!(2-TO')! ^ ^ 

('2') — 1 

The constant coefficients of the inverse matrix A^!^^^, can be derived using the inverse Lorentz transformation 
or the identity 

2 

m=-2 

For I > 2 the transformation of the harmonics is more involved. Eq. (|3.36p can easily be generalized only in the case of 
pure rotations, due to the decoupling of harmonics with different values of I, but the general Lorentz transformations 
are more difficult to work out explicitly. 

Below we discuss two special cases in which the formalism simplifies considerably: pure rotations and non-relativistic 
boosts with \V\ ^ 1. 
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1. Transformations under rotations 



Since spin-s SWSHs with a given value of I are a representation of the rotation group, for pure rotations we get 



where D^,^^ has been defined in Eq. (|3.26p . Since if = 1 for pure rotations, Eq. p.35p gives 

p2ix4(2) 



mm 



(3.40) 



(3.41) 



with ^^^)„j a constant matrix. As we discuss below, this property holds for all Vs: D^^,^ — e^^^A^^,^, and therefore 



2-' ;m' 



(3.42) 



with A\^,^ a constant matrix. 

By appropriately choosing the primed frame, a general rotation can be reduced to a rotation mapping any vector 
(such as the total angular momentum) to the z'-axis: 



(Ja;, Jy, Jz) = (Jsin6'o cos 00, •^sin6'o sin0o, '^cos^o) (0, 0, J) . 

As discussed in Appendix [b1 the stereographic representation of such a rotation has the form 



a — cqC 2 



b — — spe 2 ^ c — sqc 2 , d = cqb 2 



where we have used the short-hand notation 

Co = cos(6'o/2) , So = sin(6lo/2) . 

(2) 

The corresponding transformation matrix for / — 2 SWSHs with s — 2, A^^!^^ is 

/ 

Kpo /I (2) _ 

V 



(3.43) 
(3.44) 
(3.45) 



Co 


-2ci^so 




-2cosl 


4 \ 

-2cosi^ 






VQsqCo{1 - 2cl) 


Aslicl - 1/4) 




-%/6soCo(l - 2cg) 


1 - 6clsl 


V6soCo(l - 2co 


) 


2cosj] 


Hic'o - 1/4) 


-VQsqCo{1 - 2cg) 




-2sqCq^ 


*0 


2coSo 




2co^so 


Co' J 



(3.46) 



Here and below, whenever we write the transformation in matrix form, different rows correspond to increasing values 
of m'(= —I, . . . , Z), and different columns correspond to increasing values of m{— —I, . . . ,1). 

As a useful special case, consider the transformation of the I = 2 harmonics under a counterclockwise rotation by 
7r/2 around the a;2-axis. In stereographic coordinates this transformation reduces to 



n/2 ^ ^/2 
^/2 ^ V2 



SO Eq. (|3.34p becomes 



-(-i)"'(C- i)^+"'(C + 1)^""' 



"2-m 



(3.47) 



(3.48) 



(2) 

The corresponding coefficients A^rn'm 



,(2) 

m' m 



( I -2 \/6 -2 1 \ 

2 -2 2 -2 

\/6 -2 \/6 

2 2 -2 -2 

V 1 2 \/6 2 1 y 



(3.49) 
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The transformation under rotations of spin-2 SWSHs with Z > 2 is shgthly more involved. Under a rotation, the 
functions Z^^^^ defined in Eq. (|3.29p transform to 



r7f n 



(i + CC> 



■(aC- 



\l — 77i —n' 



[bC + d) 



-\-n 



{< + c) 



1-2-71 



{hC + df =e''^Y.P'^^mZ?ra. (3-50) 



which defines (impHcitly) the constant coefficients F^'m- l|3-50P the indices n,n' are defined in the ranges 



max(0, —2 — m) < n < I — max(m, 2) , max(0, —2 — m') < n' < I — max(m', 2) , 



(3.51) 



respectively. 

Correspondingly, the spherical harmonics transform to 



2YL' = (-1)' 



2/ + 1 {I + m'y.{l ~ m'y. 

{1 + 2)1(1-2)1 



2/ + 1 /(/ + m)!(/-TO)! 
"l^V (1 + 2)1(1 ~2)\ 



1-2 
n' 



1-2 
n 



where we defined 



A 



(0 



^ (I + m'y.(l - m'y. \ n 



1-2 



1 + 2 
I — m! — n' 

1 + 2 
I — m — n 



1 + 2 

I — m' — 



p2ix \^ pn'n yn 
/ J ni'ni hn 



m'm Im ' 



(/ + m)!(/-m)! fl-2 
n 



1 + 2 

I — m — n 



pn n 
m'm ' 



(3.52) 



(3.53) 



Remarkably, the right-hand side of Eq. (|3.53p does not seem to depend on the index n. This property implies 
Eq. (|3.42p . which shows explicitly that SWSHs with fixed I and s form an irreducible representation of the rotation 
group, as expected from group-theoretical considerations. We do not have an analytic proof that a[^^,^ is independent 
of n. However we checked this property for I < 10 and all values of m,m' (which include all cases of interest for 
gravitational wave detection), and we expect it to hold in general. 

The coefficients A^^,^ can explicitly be computed case-by-case from Eq. (|3.53p . Note that the phase factor e^'-^ is 
independent oil. As we will show, this phase factor cancels out in the transformation of the Weyl scalars. 

For reference, let us evaluate explicitly the components of the rotation matrix A'f^^,^ for I = 3. For a general rotation 
by angles (^o , M the transformation matrix A)^',^ is 



im<pQ A(-i) _ 
^ ^m'm ^ 



.0 

'IScqSq 
2V54sl 



-VGcIso ViE^sl -2V5sl4 

_c^(l-6s2) -^5/2clso(2-6sl) VSOslcKl - 2sl) 

5/_24so(2 - 6sl) 4(1 - lOsl + 15s^) - [2coSo -(- 5 sin(36'o)] /16 

30 clsl(l - 2sl) 2VtcoSo(l - 5sg -I- 5sf^) [3 cos 6*0 -f 5 cos(36io)] /8 

572cosg(4- 6sg) s^(6- 20sg + 15^) \/3 [sin^o + 5sin(36'o)] /16 

st(5 - 6sl) v^cos3(4 - 6sl) V^Oslcli^ - 2sl) 

VGcos^o VlEc^st 2V5sl4 



(3.54) 

where again we have used the short-hand notation defined in Eq. (|3.45p . The missing components can be obtained 
with the help of the general relation 



^m'm — 



m-\-m 



A 



(0 

-m' —m ' 



(3.55) 



In the special case of a 7r/2 counterclockwise rotation around the a;2-axis, from these relations and Eq. (|3.47p we get 



m'm 



















\ 




1 


-V6 


VT5 


-2%/5 


VT5 


-V6 


1 






V6 


-4 


VTo 





-VTo 


4 




1 






-Vio 


-1 


2%/3 


-1 


-VTo 


VT5 


8 




2^/5 





-2%/3 





2^/3 





-2%/5 








Vio 


-1 


-2V3 


-1 


VTo 


VT5 








4 


\/To 







-4 






I 


1 


V6 




2V5 


VT5 


^/6 


1 / 



(3.56) 
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2. Transformations under non-relativistic boosts 

In the case of boosts, Appendix iBl shows that the conformal factor is K = 1 + Vifii to hnear order in the boost 
velocity. At the same order, the transformation in the stereographic representation is described by the parameters 

a-l-y, '^^^ + Y' (^-^^^ 

where we have defined V± = Vi :k iV2, and 1^ <C 1 is the boost velocity along the Xi-axis (see Appendix [B|) . The 
components A^J^ can now be determined from Eqs. (|3.34p and (|3.37p . To first order in Vi they read 



4(2) 

rn'm 



/I 


-2^3 


v+ 








\ 




V- 




y/3/2V+ 















1 


\/3/2V+ 















I + V3 


v+ 


V 











V- 


1 + 2V3J 



(3.58) 



We stress that, although only I — 2 harmonics formally appear in the expansion (j3.35p . for a boost there is no true 
decoupling in I, since the conformal factor K has a non-trivial dependence on 9 and 4>. Harmonics with I > 2 would 
appear if we expanded in constant coefficients, as in Eq. (|3.22p . 

C. Transformations of the Weyl scalar 

In order to find the Lorentz transformation of the Weyl scalar ^'4, as given in Eqs. (|3.23p and (|3.24p . we need to 
determine the inverse transformation of spin-2 SWSHs, i.e. the coefficients D^rn i'm'- ^^^^ ^° ^^'^ following. 

1. Transformations under rotations 
Under rotations, using Eqs. (|3.25aP and (|3.42p . Eq. (|3.23p reduces to 

^I'iu', 9', 0') = e-'^" 2^-^' 211™' {e\ , (3.59a) 

Im' 

2^L' = E 2^-;™^!!"' • (3.59b) 

m 

It is worth noting that the phase factors e^^'-^ cancel, and the transformation under rotations of the coefficients of 
the multipolar expansion of ^4 coincides with the inverse transformation of spin-2 SWSHs. 

The transformation laws for the standard decomposition of the Weyl scalar ^1^4 in s = — 2 harmonics can be obtained 
with the help of relations (|3.5p . (|3.7p . However, a simpler relation emerges if one uses Eq. p.55p . Using the complex 
conjugate of Eqs. (|3.59bp and p.55p one can show that 

^l'{u', 9', 0') - e''^" E -2^-^' -2Yim' {9', <P') , (3.60a) 

Im' 

-2i^'lm' - E -2'/';mA!iL"' • (3.60b) 
m 

This means that the coefficients of s = — 2 and s = +2 transform with the same representation. Some applications of 
Eq. (|3.60bp to numerical relativity simulations of black hole mergers are discussed in Section IIVI 

2. Transformations under boosts 

The case of pure boosts is more complex than a pure rotation. The conformal factor K is now a function of 9 and 
0, so it is not trivial to extract the transformation law of the coefficients 2'0;Tn from Eq. ()3.23p . 

Quite interestingly, the factor K{9, (j)) couples harmonics with different values of /. To see how this coupling comes 
about and how boosts affect the mode decomposition, let us consider a gravitational wave which in the original frame 
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consists of an (Z = 2, m = 0) mode only. If we move to a frame which is boosted by V3 along the Xs-axis and we 
neglect the phase factor e~'"R'" (e.g. by considering wavefronts for which u' = 0), from Eq. (|3.23p we get 



*°'(?/', e', « (1 - 51^3 COS0') 2^20 4m' ' ^¥2^' (9', <f>') 



(3.61) 



where we expanded K ^ = (1 + V3 cosO) ^ ~ 1 — 5V3 cos^. The inverse transformation matrix, ignoring second-order 
terms, is 



,(2)-i 



/ 1 - 2V3 

1 + V"3 

10 

I-V3 



V 











(3.62) 



1 - 2^3 / 



If we equate Eq. ([^31]) to Eq. (|3.59ap we get 



E 



2V^L 2>lm(f?', 0') = 2^20(1 - 5^3 COS^') 2^20(0', (j)') = 2V^20 



ji"2o(e',0') 



V7 



1^3 2>30(e'>') 



where we used the property 



cose' oY-i 



2-120^ 



%/7 



2^30 ( 



(3.63) 



(3.64) 



Therefore, the boost produces an I — i component: the nonvanishing components of the multipolar expansion are 

2V'20 = 2V'20 , 



2^30 = 7E^3 2^^: 



20 • 



(3.65) 



An I = 2 waveform is seen as distorted in the boosted frame, with contributions also from I = i. If we took into 
account higher powers of V we would find that higher multipoles would also be excited. 

Some considerations on the physical relevance of this coupling are in order. For unequal-mass spinning black hole 
binaries, the final black hole can be kicked with velocities as large as 4000 km/s [l^, 0, [1^, [l^l- In the center-of- 
mass frame, numerical simulations (and possibly also post-Newtonian expansions) will see an ? = 3 component due 
entirely to the fact that the black hole is boosted. The I — 3 component can be as large as V3 ^ 10^^ times the 
I = 2 component. Boost-induced multipole coupling will certainly be relevant in small impact-parameter collisions of 
ultrarelativistic black holes: in this case the velocity is a sizeable fraction of the speed of light, and the energy transfer 
between modes can have important consequences. A detailed study of this problem will be reported elsewhere. 



IV. APPLICATION TO NUMERICAL SIMULATIONS AND BLACK HOLE PERTURBATION THEORY 



In this section we apply the tools developed in the rest of the paper to numerical simulations of binary black-hole 
mergers and to a classical problem in black-hole perturbation theory. We examine in detail two scenarios. 

In the first scenario, we consider two possible physical situations: the head-on collision of equal-mass, non-spinning 
black holes in numerical relativity, and the infall of a point particle into a Schwarzschild black hole in perturbation 
theory. We study the collision in two different frames: (1) a frame where the collision (or infall) occurs along the 
X-axis, making maximal use of the symmetries of the system, and (2) a frame where the collision occurs in the 
equatorial plane (along the x-axis, for the equal-mass case). The two frames are obviously related by a 7r/2 rotation: 
^ 9' ^ 9 — IT /2. Our main purpose is to show in a relatively simple but non-trivial setting that our formalism is 
consistent with the predictions of numerical simulations and of black-hole perturbation theory. 

In the second scenario we consider a more general case: the merger of unequal-mass, spinning black holes in a so- 
called "superkick" configuration [H, [23, H^, Hj] ■ This configuration is general enough to provide a good case study for 
generic, unequal-mass mergers of spinning black holes. Our main goal in this exercise is to stress the importance, in the 
context of generic mergers, of relating two reference frames: the "simulation frame" (where the binary initially orbits 
on the x — y plane, so that the orbital angular momentum initially points in the z-direction) and the "remnant frame" 
(defined as the frame where the final black hole's spin is aligned with the z-axis). To a good level of approximation, 
if we ignore the (relatively small) boost, the two frames are related by a rotation. We suggest a simple method to 
estimate the angles (^0 i 4>q) describing this rotation. Then we argue that the analysis of merger waveforms is more 
"natural" in the remnant frame. More specifically, we show that (i) the merger radiation, as seen from the z-axis, 
is circularly polarized in the remnant frame, but not in the simulation frame, and (ii) ringdown frequencies provide 
reasonable estimates of the final black-hole spin in the remnant frame, but not in the simulation frame. 
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A. Head-on collisions of non-spinning black holes 

Let us start by discussing how rotations affect the description of two "simple," highly symmetric processes leading 
to the generation of gravitational waves: the head-on collisions of equal-mass, non-spinning black holes in numerical 
relativity, and the infall of a point particle into a Schwarzschild black hole in perturbation theory. 

1. Equal mass head-on collisions in numerical relativity 

For an equal mass head-on collision of non-spinning holes along the z-axis, by symmetry, only modes with m ~ 
are excited (see [i^ and references therein). In Ref. [4^ we studied the transition from circular inspirals to head-on 
collisions by considering a sequence of eccentric inspirals. Initially the black holes were placed on the x-axis. Along 
the sequence, we fixed the system's binding energy. The (equal and opposite) linear momenta on the two holes 
were directed along the y-axis, and gradually reduced from the value P = Pcirc required to produce a quasi-circular 
inspiral down to P = 0. In the limit P — the black holes collide head-on along the a;-axis. By analyzing the 
resulting radiation we observed that it contains components with (/ = 2,m = ±2) and {I = 2,m = 0), and that 
these components have roughly comparable amplitude (see top left panel in Fig. 2 of ^46] ) . Now we are in a position 
to predict quantitatively the relative amplitude of these components. In the frame where the holes collide along the 
2-axis, we should only have a component with (I = 2,m = 0). By colliding the holes along the x-axis, we are effectively 
performing a 7r/2 rotation in the ^-direction. But then from Eq. (|3.49p we immediately see that the (/ = 2,to = 2) 
and (Z = 2, TO = 0) wave amplitudes should be proportional to each other, with a proportionality factor of — \/6/2. In 
Fig. [T]we show that this prediction is in excellent agreement with the numerical simulations. The error (short-dashed, 
green curve) is orders of magnitude below the wave signal except for those times when the GWs are so weak that ^'4 
is dominated by numerical noise; that is, during the early inspiral, once the junk radiation has left the system, and in 
the late ringdown stages. Ref. ^] independently derived a similar result in the context of head-on collisions of boson 
stars. 




t/M 

FIG. 1: Check of the transformation law for the I = 2 modes of the head-on collisions studied in [4^. To better quantify the 
error, the short-dashed (green) line shows |Re[M(_2V'22 -I- V6_2''/'2o/2)] | . This quantity should be zero, and indeed it is typically 
much smaller than the corresponding amplitudes, with the exception of the early inspiral and late ringdown phases. 



2. Infall of particles in black-hole perturbation theory 

The previous discussion referred to equal-mass black hole collisions. The limit where one black hole is much smaller 
than the other (Mi ^ M2) is also of great interest, since it can be described within black hole perturbation theory. 
In this framework the small black hole is considered as a perturbation in the geometry of the large black hole, and 
one can expand Einstein's equations to first order in the parameter M1/M2. The resulting system of equations can 
be further reduced to a single ordinary differential equation with a source term C The functional form of the source 
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term depends of course on the frame chosen to do the calculation. In particularly symmetric situations (e.g. for the 
radial infall of a small body into a more massive black hole) the source term has a simple analytical form. For this 
reason, the radial infall of a point particle into a black hole is a good testing ground for the formalism we developed in 
this paper. We consider the source term in two frames: one where the infall occurs along the symmetry axis (so that 
only modes with m — are non- vanishing) and a frame rotated by 7r/2, so that the collision occurs in the equatorial 
plane. 

Details of the formalism are discussed in Appendix|Dj Here we only report the main result: for particles falling along 
the axis and in the equatorial plane the source terms are given by Eqs. (|D10P - (jDlip and (jD14aP - (jD14bp . respectively. 
Since we work in linear perturbation theory the waveforms are proportional to C, and therefore we infer that 

.equator /F „;,axis .equator rr" ,/,axis a 

/equator o ' /equator ' /equator A/ o ' /equator /p ' V ■ / 

-2'lp20 ^ -2^20 -2V'3±1 ^ -2V'33 

These proportionality relations are perfectly consistent with the transformation matrices p.49p and (|3.56p . 



B. Generic rotation angles: unequal-mass "superkick" runs 

As a non-trivial practical application of the transformation laws (j3.46p and (j3.60bp . in this section we study the 
inspiral and merger of unequal-mass, initially spinning black holes. Our study is a rather striking demonstration of 
the importance of choosing the "right" reference frame to analyze the radiation produced by a merger event. 

Our case-study is the "sk4" run described in Ref. [l^. Our black- hole binary has mass ratio q = Mi /A/2 — 4. The 
Kerr parameters of both holes have the same magnitude |ji| = |j2| = 0.727, but opposite directions, and lie on the 
orbital plane at the beginning of the simulation: this is sometimes referred to as a "superkick" configuration. 

Consider the dimensionless spin vector of the remnant black hole J — {Jx,Jy,Jz)- Define Jxy — (Jx + CfyY^"^ ^ 
and let J = {J^ -\- Jy -I- J'^Y^'^ be the modulus of the final spin. Simple flat-space geometry suggests that a rough 
but reasonable estimate of the rotation angles relating the simulation frame to the remnant frame should be 

6*0 = arcsin(J'^y/J') , = arccos ( J'^,/ J'^jj,) . (4.2) 

For our sk4 run, the final spin can be obtained from angular momentum balance (see Appendix[C]for details). The final 
spin magnitude estimated in this way \s J = 0.642, and the individual spin components are Jx = 0.424, Jy — 0.036, 
Jz — 0.480. The system experiences a significant realignment: the estimated rotation angles are — 41.5° and 
(/)o = 4.8°. Now we can use this estimate of the rotation angles, together with relations p.46p and (|3.60bp . to 
transform the "simulation frame" waveforms to the "remnant frame" waveforms. 




{ty/M (ty/M 

FIG. 2; Top: Quadrupolar {I — 2) component of |Re(Mr^'4)| and |Im(Mr 'I'4)| as seen by an observer located on the z-axis, 
for the sk4 run of Ref. 19]. Bottom: degree of elliptic polarization of the waveforms (compare Fig. 28 in Ref. [T^). The left 
panel refers to the simulation frame, the right panel to the remnant frame. 

Our results are summarized in Fig. [21 The upper left panel of Fig. [2] shows the dominant, quadrupolar {I = 2) 
component of ^'4, as seen by an observer located on the z-axis in the simulation frame. To obtain this quantity we 
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take the sum 

^ -2V'i™(i,r)_2F/™(0,</'), (4.3) 

1=2 m=-2,0,2 

of alH = 2 components with different values of m (i.e. m = 2, m = and m = —2), and we evaluate it for 9 = (j) = 0. 
The upper right panel plots the quadrupolar (I = 2) component of ^'4, as seen by an observer on the z-axis in the 
remnant frame. Again, this quantity is obtained by taking the sum of all Z = 2 components in the new frame, as 
prescribed by the transformation laws for 1 = 2. In both plots the time axis is shifted by Iq, the time corresponding 
to the peak in the waveform's amplitude, which is also a rough measure for the formation of a common apparent 
horizon and the beginning of ringdown (see e.g. [3]). 

Some simplification in the waveform is clearly seen in the remnant frame, especially in the late (ringdown) phase. 
In the remnant frame, the ringdown part of the radiation is well described by a single, "clean" damped exponential, 
and the envelope of the maxima is well fitted by a straight line. In the simulation frame the envelope of the maxima 
is an oscillating function. This happens because the ringdown radiation is effectively a sum of components with the 
same l(= 2) but different values of m, and this produces beatings between different quasinormal modes. 

An alternative measure of the improvement produced by using the remnant frame is given by the degree of elliptic 
polarization Pe, shown in the bottom panels of Fig. [51 This quantity was defined and discussed in Appendix D 
of Ref. [1^. It has the property that Pe = 1 for circularly polarized waves, and Pe = for linear polarization. 
Fig. [2] shows that the ringdown signal is not circularly polarized as seen from the z-axis in the simulation frame, but 
it is indeed circularly polarized (to a very good approximation) in the remnant frame. This is a rather impressive 
demonstration that the merger-ringdown waveforms look simpler in the remnant frame. On physical grounds this 
makes perfect sense, since the remnant frame (unlike the simulation frame) is well adapted to the final hole. For 
generic precessing binaries, it is still possible that switching to the remnant frame does not simplify the early inspiral 
waveform. A detailed analysis of longer, generic inspiral simulations will be a topic for future research. 

It was observed a lon g ti me ago that the polarization state of the waves is a useful discriminant between alternative 
theories of gravity [11,14^. More recently, it has been suggested that LISA could place interesting bounds on the 
amplitude birefringence (i.e. the polarization dependent amplification, or attenuation, of the wave amplitude with 
distance) predicted by Chern-Simons gravity [5Q] . The fact that the merger-ringdown waveform is circularly polarized 
in the remnant frame may have useful applications to test general relativity using strong-field observations of merging 
binaries. 
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FIG. 3: Final angular momentum from the {I = 2 ,m = 2) mode in the remnant frame (solid line) and in the simulation frame 
(dotted line). Dashed thick (thin) red lines: we arbitrarily increase (decrease) 4>o by A(^o = 4°. Dash-dash-dotted thick (thin) 
blue lines: we arbitrarily increase (decrease) 9o by A^o = 8°. The horizontal line represents the estimate of the final spin 
[J' = 0.642) obtained from balance arguments. 

As an additional demonstration of the usefulness of the remnant frame it is useful to look at the ringdown radiation. 
We use a Prony analysis of the I — m ^ 2 ringdown mode to extract the quality factor of the quasinormal ringing 
[33 |. From the quality factor we compute the Kerr parameter of the final hole (see [16] and [5l| for further details), 
comparing this estimate with the prediction from angular momentum balance arguments. Our results are shown in 
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Fig. [3l where we show the spin estimated by considering different portions of the ringdown signal, starting at time 
t > Iq. In the remnant frame, the final spin estimated from the ringdown signal is much closer to the value estimated 
from balance arguments. The waveform in the simulation frame contains a mixture of various modes of comparable 
amplitude and different damping times, and this produces a large error in the estimation of the final spin. 

As an exercise, in Fig. [3] we also show the final spin estimated in different frames which are "reasonably close" to 
the remnant frame, considering frames that deviate from our "best guess" for the simulation frame by small angles. 
For illustration we plot spin estimates obtained assuming Oq = 41.5° ± A6'o or (pQ = 4.8° ± A(/)o, where A0o = 4° 
and A6'o = 8°. This shows that the good results we obtained in the remnant frame are not accidental: for example, 
adding (or subtracting) A6'o = 8° to our "best estimate" for the rotation angle 6*0 seems to produce a systematic 
underestimate (or overestimate) of the final spin, as obtained from angular momentum balance arguments. Systematic 
errors in our numerical simulations are admittedly quite large, but the significant bias observed for "incorrect" values 
of A^o is an indication that frames deviating from the remnant frame by small amounts of rotation can yield poor 
agreement with the final spin estimated from balance arguments. 



We have provided explicit, ready-to- use formulae to transform multipolar components of gravitational radiation 
under rotations and non-relativistic boosts. For generic spinning binaries our results indicate that gravitational 
waveforms are simplest in the "remnant frame" (see Figs. [5] and [3]). We provide a practical recipe to identify this 
frame, and we suggest that all results from multipolar analyses of merger simulations should be expressed in the 
remnant frame to simplify comparisons. Our work will be useful to standardize file formats (soj and for the injection 
of numerical waveforms in the data analysis of gravitational- wave detectors [3l| . Other interesting applications will be 

(1) an analysis of the multipolar coupling induced by large boosts in the collision of ultrarelativistic black holes, and 

(2) a study of the possibility to determine the spin orientation of a binary merger remnant using gravitational-wave 
observations on Earth and in space. Both topics are currently under investigation. 
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Transformations of the form (|2.2p and (|2.4[) are conformal transformations of the (6', (/))-sphere into itself, with 
conformal factor K . They correspond to Lorentz transformations of spacetime; indeed, the conformal group SL{2, C) 
and the Lorentz group S0{3, 1) are isomorphic, and any Lorentz transformation on spacetime has, on the sphere, the 



In this Appendix we show this explicitly for the case of a non-relativistic boost with three-velocity Vi such that 
\V\ ^ 1. We also assume that, since we are interested in asymptotic transformations, \u\/r ^ 1. We introduce 
polar coordinates such that Xi — rrii — r(sin6' cos sinf? sin0, cos6'), and neglect terms of order 0{V'^) and of order 
0{\u\/r). Then under a boost 
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APPENDIX A: BOOSTS AS BMS TRANSFORMATIONS 





t - ViXi . 



(Al) 
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In polar coordinates we have 



It follows that, to linear order in V, 



t' = t-rViU,. (A2) 



u' ^t' -r' = u{l + V,n,) = Ku , (A3) 
with K{9, (()) — 1 + Villi. This is our main result. Neglecting terms of order 0{\u\/r), 

r' = r{l- !ii^K;n»^ = K'^r . (A4) 

The transformation of the director cosines, i.e. of {9, (f>), can be found as follows: 

x[ = r'n^ — (r — tVjnj)n[ — Xi — Vit — rrii — Vit , (A5) 



therefore 



/ ^ \ ^ / U ^ T \ 11 T 

n\=[l + -V,n,j - -V, = [l + ^V,n, j - = (1 + V,n^ 



n,-V,^ Kn, - V, . (A6) 



These are well-defined director cosines: indeed, to linear order in V, 

n'X - [Kn, - Vi){Kn, - Vi) = - 2n,V, = 1 . (A7) 

Using the equality ^Kn^Vi = 2(1 + VinijriiVi = 2niVi + 0{V'^) and consistently neglecting terms of 0{V^), it can be 
shown that the metric on the two-sphere changes as follows: 

d9'^ + sin^ 9'd(t)'^ = dn'^dn'i = K'^dn.dm = {d9^ + sin^ 9d(j)^) . (AS) 

The boosts, then, act on {u,9,(j)) as in Eqs. (|2.2p and (|2.4p . Notice that the flat spacetime metric does not change 
under a boost, i.e. 

ds^ = —du^ — dudr + r^duidfii — —du'^ — 2dudr' + r'^dn'idn'i , (A9) 



but this is not obvious by replacing (|A3p , (jA4p and (jA6p in the above relation. The reason is that we have neglected 
terms with u/r <C 1. As a qualitative check of the expression K ~ I + ViUi, let us consider the case in which 
the observers are boosted towards the source. Then ViUi < 0, and K — 1 + Viiii < 1. Therefore, the observed 
gravitational-wave frequency changes as uj' — to/ K > uj, i.e. we see a blueshifted wave, as expected. 

We have shown that a boost with velocity V such that |V^| <C 1 has conformal factor K — 1 + ViUi. It affects the 
transformation of the asymptotic Weyl scalar ^'4, Eq. (|3.21|) . by terms linear in V. We also have to consider the effect 
of the change in the coordinates, which is of order V as can be seen from (|A3|) . (jA4l) and (|A6|) . and is harder to derive 
explicitly. It is derived in a particular case, in Section IIII C 21 Notice that, if the black hole kick determines both 
a boost and a translation, the translation only appears in the change of coordinates, and not in the transformation 
(I3.2ip . A translation is expected to give a contribution of order |(52;|/r, which is extremely small far away from the 
source: for the purpose of detection, the boost is the most relevant transformation. 

APPENDIX B: THE STEREOGRAPHIC REPRESENTATION 

Conformal transformations of the sphere (which correspond to asymptotic Lorentz transformations of spacetime) 
can be expressed as conformal transformations on the complex plane obtained by projecting the sphere from the north 
pole to the equatorial plane (2/1,2/2)- From the geometry in Fig. [H defining a complex variable ( — yi + iy2 we have 

C = cot^e'^ (Bl) 

The 2-sphere metric can be written as 

d9'^ + sin^ 9d(l>^ = 4(1 -I- CCy^dCdC . (B2) 
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The director cosines — (sin 6* cos 0, sin 6* sinc/i, cos 0) are related to C by 



2C CC-1 

n\ + 1719 = , n-x = , \Do) 

CC + 1 CC + 1 



and by the inverse relation, ^ = (rii + in2)/(l — f^s)- An SL{2, C) transformation 

maps C to C as in Eq. (|2.6p . with ad — be — 1. Notice that if ad — foe 7^ we can rescale the four elements without 
changing the transformation, in order to have ad ^ be = 1. This can be shown to be a conformal transformation, i.e. 
the 2-sphere metric changes as in Eq. (j2.2p . which in stereographic coordinates means 

4(1 + CCT^dC'dC' - 4(1 + CO^'^^C^^C ■ (B5) 

A straightforward calculation shows that 

1 ^ jaC + c){aC + c) + {bC + d){bC + d) 

CC + 1 ■ ^ ' 

The phase x appearing in Eqs. (|3.15p and (|3.17p is given by (see Eq. (3.12) in 0) 



dC'/dC [bC + df 



e-2ix^^w^= v^^^:^. (B7) 



1. Rotations in stereograpiiic coordinates 

The main purpose of this Section is to derive the specialization of ()2.6|) to a rotation by {9q,(I)o), i-e. Eq. (|3.44|) . 
and its special case Eq. (|3.47p for rotations by an angle 6*0 = tt/2. Rotations are SL{2, C) transformations for which 
K = 1. Using Eq. (|B6p this implies 



(aC + c)(aC + c) + (6C + d)(bC + J) = CC + 1 , (B8) 

and therefore 

\a\' + \b\' = 1, 
|cp + |dp = 1, 

ac + bd = 0. (B9) 

But ad - be ^ 1, so A e SU{2) and 

a = d, b=-c. (BIO) 

Let us first consider, for simplicity, a rotation of tt/2 around the X2-axis; it maps X2 X2, xi — > X3 and X3 —xi, 
i.e. 

sin 9' cos (j)' — — cos 9 , 
sin 0' sin (j>' = sin 9 sin (f) , 

cos 6*' = sin 6* cos 0. (Bll) 
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In particular, the plane </) = </)' = is left unchanged, and 0' = 9 — ^. The corresponding function C'(C) niaps (in the 
plane = 0) —1 in 0, in 1 and 1 in oo, therefore it must have the form 

^, _ aC + c _ ajC + l) , . 

^ -6CTd-a(~C+l)' ^ ^ 

and imposing ad — be = 1 this gives Eq. (|3.47[) . As expected, this transformation satisfies the relations (jB9p . thus 
K = 1. 

By appropriately choosing the primed frame, a general rotation maps a vector with generic orientation to the z'-axis: 
(•/x, ^y, ^z) = (>/ sin 6*0 cos 00, J sin 6*0 sin 00, J cos ^o) — > (0,0, J). (B13) 
Therefore, it maps Co = cot(6'o/2)e''^»/2 ^' ^ Imposing (|B9|) (i.e. K ^ 1) yields 

-sin|ie-'^ +cos%e'T^ 

i.e. Eq. ^ZM- 

2. Boosts in stereographic coordinates 

Let us now express boosts in stereographic coordinates. As shown in Appendix El the conformal factor of a boost 
with velocity Vi such that |y| ^ 1 is if = 1 + T^n.i. We define 

V±^Vi± iV2 , (B15) 

thus 

V.n. = V^"^^^ + vil^l±^ + ^^3^3 = ^ + V^-C + ^3(CC - 1)] . (B16) 
Therefore, to linear order in V ^ 

K = l + ^^[V+C + V-C + V^{a-l)\ , (B17) 

K'' - 1-^^[V^+C + T^-C + ^3(CC-1)] . (B18) 

Now we will work out the expression of a boost as an SL{2, C) transformation. A four-vector in spacetime can be 
written as an SL(2, C) matrix X: 

X = a,x,= f ^"+."^ -^~'^-A (B19) 

where cro is the identity matrix and <Ti {i = 1,2,3) are the ordinary Pauli matrices. A boost in the xi direction is 
represented by the matrix 

io u ( cosh^ -sinh^\ 

^ = --a, smh - = _ ^^^^ J j , (B20) 

where the velocity of the boost is y = tanhw. Indeed, the transformed vector is 

^^^1 _ f xo coshw — xi sinhw + — xq sinh w + xi coshw — ix2 \ ('B21) 
I — xq sinhw + xi coshw + ix2 xo coshcj — xi sinhw — X3 J ' 

where a dagger denotes Hermitian conjugation. For a boost in a generic direction, defining three parameters uii 
[i = 1, 2, 3) we have 

A = e ^ = cosh - f i sinh ^ . (B22) 
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To linear order lu = V + ©(V^), and 



2 ' 2 

which is the result quoted in Eq. (|3.57p . Such a boost maps ( — cot |e"^ into 



C' = ^= \ . ,r v_ 1-^ + ^ (B24) 



bC + d i + 

and 

+ 2 2 

As a check of Eq. (|3.57p it is easy to verify that the inverse conformal factor (jB6|) , when expanded to linear order in 
V, agrees with Eq. (|Bl8l) . 



APPENDIX C: 



CALCULATION OF THE RADIATED ANGULAR MOMENTUM IN NUMERICAL 
EVOLUTIONS OF BLACK-HOLE BINARIES 



All numerical black- hole binary evolutions in this paper have been performed using the Lean code 43]. We 
compute the energy, linear and angular momentum radiated during merger using the Newman-Penrose scalar 5*4. 
The calculation of ^'4, as implemented in the Lean code, is described in detail in Appendix C of Ref. 'A3\. Here we 
assume 'if 4 as given on an extraction sphere of constant coordinate radius r^^. 

The radiated energy, the {x ,y , z) components of the radiated linear momentum and the z component of the radiated 
angular momentum are then given, for example, by Eqs. (22)- (24) of Ref. [52|]. For the simplest binary black hole 
configurations these quantities encapsulate the entire information about the radiated momenta, because the x and y 
components of Jrad are known to vanish from symmetry arguments. However, for more general classes of spinning 
binaries (such as those considered in this paper) we can no longer neglect these components. The calculation of J^j^rad 
and Jy^rad from ^'4 has recently been discussed in [sl,!!^!- Here we follow Eqs. (2.24, 2.25) of [13] and calculate the 
radiated angular momentum from 



dJ i I- 



ad 



dt 



where 



lim Re (p Ji 




-^Adt dt 



-^Adt dn 



(CI) 



sin dJx = — sin 9 sin (f)dg — cos 9 cos 
sin 9Jy = sin 9 cos (pdg — cos 9 sin 



+ IS cost 
is cos 6, 



(C2) 



and the spin weight s — ~2. 

As a test of the numerical implementation of these expressions we followed the suggestion of Sec. HI of Ref. [s^ . 
We evolved the binary configuration labelled "Rl" in table I of Ref. [55;] three times, the orbital angular momentum 
L being aligned with the z, y and x axis, but always using the z axis as reference for the wave extraction. We use 
the grid setup of Table I in [43] adapted to the orientation of L, using a comparatively low resolution h — 1/32 for 
these test simulations. The resulting radiated angular momentum is shown in Fig. [51 The radiated angular momenta 
in the three simulations agree within a few percent, which for these low resolution runs is less than the uncertainties 
resulting from the discretization. 



APPENDIX D: RADIAL INFALL OF PARTICLES INTO NON-ROTATING BLACK HOLES 



Small perturbations of the Kerr geometry are usually described using the Teukolsky formalism fJT, '56*] , but some 
practical applications require a modification of the equations. The Sasaki-Nakamura formalism ,57i . i58. . i59l . [60l . [6ll . 
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FIG. 5: Radiated angular momentum as a function of time from the Rl simulation with the orbital angular momentum aligned 
with the X, y and z axis (see text). 



|62| . l63l |6J, |65| is one such modification that is very useful (among other applications) to deal with the infall of 
point particles into Kerr black holes. For simplicity, here we consider a highly relativistic point particle falling into a 
non-rotating black hole, since in this case the equations simplify considerably (66l . [67j . Our purpose is to check the 
formalism developed in the text by considering the same problem in two different reference frames. In the first case, we 
consider the radial infall of a particle along the z-axis, making full use of the symmetries of the problem. In the second 
case, we deliberately let the particle fall in along a radial trajectory on the equatorial plane, destroying the symmetry 
of the problem. The point of this Appendix is to derive the source terms of the Sasaki-Nakamura equation in these 
two configurations; using such expression, in Section llV A 21 we show that the transformation properties derived in 
this paper are perfectly compatible with the transformation laws of the source term of the Sasaki-Nakamura equation. 

Details about the Teukolsky formalism may be found in the original literature [4l|, [s^ , and also in [6^ . Working 
in the Kinnersley tetrad, the equations for the Newman-Penrose quantities decouple and separate, giving rise to the 
Teukolsky equation for a radial function R [4ll . [56j : 

A-'^(a'+^^r)- sVR=-T, (Dl) 
ar \ ar J 

where we defined A = — 2Mr , and M is the mass of the black hole. The quantity T is a source term due to the 
stress-energy of the particle, and the potential 

-K^ isKA' 2isK' , , , , 

^V^^ + — —-2isK'+A. (D2) 

We work in the frequency domain, assuming a time dependence of the form e^"^*. The function K = r'^uj, and primes 
denote derivatives with respect to r. The quantity s is the spin weight (or helicity) of the field under consideration (see 
for example [39|). and m is an azimuthal quantum number. We are interested in gravitational perturbations, which 
have spin- weight s = —2. For non-rotating black holes the separation constant = + — s{s + l). The source term 
T appearing in Eq. (|Dip depends on the stress-energy tensor of the perturbation under consideration (see e.g. [gsI]). 
For a test particle falling in from infinity with zero velocity T ~ r^/^, and the potential sV is always long-ranged. 
This implies that when we try to solve the Teukolsky equation (jPip numerically we will run into problems. Usually 
the numerical solution is accomplished by integrating the source term multiplied by certain homogeneous solutions. 
At infinity the integral will not be well-defined: the source term blows up there, and because the potential is long 
range, the homogeneous solution will also blow up at infinity. We can remedy this long-range nature of the potential, 
and at the same time regularize the source term, by using the Sasaki-Nakamura formalism. After some manipulation, 
the Teukolsky equation (jPip may be cast in the Sasaki-Nakamura form [s^, [H, |H, HO, [U, IH, [g^l : 

-^X{ij,r)-J'^X{u;,r)-UX{ij,r)^C. (D3) 
art dr^ 
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The functions T and U can be found in the original hterature [57|, l58|, l59|, l60|, l61|, [62, |63| ■ The source term L is given 

by 



(D4) 



and W satisfies 



W" = -—Texp 



(D5) 



The tortoise coordinate defined by dr^jdr — r"^ / IS., ranges from — oo at the horizon to +00 at spatial infinity. The 
Sasaki-Nakamura equation (jD3p must be solved imposing an outgoing-radiation boundary condition at infinity: 

r» — > 00. (D6) 



X(tj,r) =X°"*e''^''- , r, ->oo. 
The two independent polarization modes of the metric, and /ix, are then given by 



i/iv = 



g ^+00 ^ 



A(A + 2) - 12iw 



(D7) 



Using a formal Green's function solution it is easy to see that X°^^ (x C. 

Let us now consider the radial infall of an ultrarelativistic point particle into a non-rotating black hole. The 
dependence of the emitted gravitational wave power on the trajectory of the particle comes entirely from the source 
C in eqs. (|D3p - (|D4|) . which in turn depends on the source term T in Teukolsky's equation (|Dip through Eq. (|D5p . 
To compute the explicit dependence, one would have to consider generic geodesies in the Kerr geometry 69], which 
in general gives an end result not very amenable to analytical treatments. However in the highly relativistic regime 
the equations simplify enormously fro . [67| , and in particular for radial infall (along the symmetry axis or along the 
equator) it is possible to find a simple closed-form analytic expression for the function W in Eq. (jD4p . Accordingly, in 
the following we will provide a closed- form expressions for C in two special situations: radial infall of an ultrarelativistic 
particle along the symmetry axis, and radial infall of an ultrarelativistic particle on the equatorial plane. We will 
further specialize to non-rotating black holes, where analytic expressions for the spin-weighted harmonics are known. 



1. Particle falling along the symmetry axis 

Consider a particle falling radially along the symmetry axis of a Kerr black hole. In this case the geodesies, written 
in Boyer-Lindquist coordinates, are 

dt r2 fdrV ^ A d(l) , , 

where the parameter eq is the energy per unit rest mass of the infalling particle. On considering highly relativistic 
particles, Cq 00, the source-term function C reduces to [66l. Igt} 



C = -i^^^e-^-'. (D9) 



Here C = 8 (_2^1o/ sin 0) g^^, the function 70 = "fo{r) can be found in [57|, |58|, |59|, |6d, l61|, 162, l63| , and ^ is the mass 
of the particle. We find 

^axis ^ ^-Ar-^ e-'"'^*, {1 = 2), (DIO) 

^axis ^ V — e-'""'*, (; = 3). (Dll) 



24 



2. Particle falling on the equatorial plane 



Now suppose that the particle is faUing radially on the equatorial plane of a Kerr black hole. On considering highly 
relativistic particles we find that the source term takes again a very simple form [6^ [67| : 



(D12) 



where 



S 



^2Yira (7r/2,0) + m_2i^™(V2,0) , 



(D13) 



and a prime stands for derivative with respect to 9. For the dominant multipolar components we get 



/^equator 

■'-'20 ~ 



^equator 
•'-3±3 



±- 



^equator 

'-2±2 — 



^equator 

^3±1 — ^ 



4w^r' 
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